MKC  Technical  Summary  Report  #2236 

MOLECULAR  TRANSPORT  EFFECTS  IN 
TURBULENT  DIFFUSION  FLAMES  AT  MODERATE 
REYNOLDS  NUMBER 

R.  W.  Bilger 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 


July  1981 


(Received  February  18,  1981) 


6 


Sp^PPred  by 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Approved  for  public  release 
Distribution  unlimited 


81  9  08 


.‘Jr-.-  *-  - 


UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 

MOLECULAR  TRANSPORT  EFFECTS  IN  TURBULENT  DIFFUSION 
AT  MODERATE  REYNOLDS  NUMBER 

R.  W.  Bilger 

Technical  Summary  Report  #2236 
July  1981 

ABSTRACT 

In  flames  molecular  dif fusivities  are  enhanced  by  the  high  temperatures 
and  can  be  of  the  same  order  as  turbulent  dif fusivities  in  flames  of  moderate 
Reynolds  number*  A  perturbation  analysis  is  used  to  quantify  effects  in  a 
hydrogen/air  diffusion  flame  which  arise  from  differential  molecular 
diffusivities.  The  analysis  uses  perturbations  about  the  equal  dlffusivity, 
adiabatic,  equilibrium  theory  commonly  used  and  yields  solutions  for  the 
average  and  higher  moments  for  the  departures  in  normalized  element  mass 
fractions  and  enthalpy.  The  results  are  compared  with  the  laser-Raman 
measurements  of  Drake  et  al.  in  a  relatively  low  Reynolds  number  flame. 
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MOLECULAR  TRANSPORT  EFFECTS  IN  TURBULENT  DIFFUSION  FLAMES 


AT  MODERATE  REYNOLDS  NUMBER 
R.  W.  Bilger 

Introduction 

In  theoretical  models  [1.3]  of  turbulent  diffusion  flames,  great  applications 
are  Introduced  by  the  assumption  that  the  mass  diffusivlties  of  all  species  are 
equal  and  equal  to  the  thermal  diffusivity,  that  Is,  all  Lewis  numbers  are  unity. 

It  can  then  be  shown  [4]  that  the  elemental  composition  and  enthalpy  at  any  point  in 
the  fluid  is  the  same  as  that  obtained  by  mixing  a  mass  5  of  the  unmixed  fuel  with 
a  mass  1-5  of  the  unmixed  oxidant.  The  mixture  fraction  5  »  5  (jt.t)  is  then 
the  only  variable  needed  to  describe  the  state  of  mixing.  No  differential  (or 
"preferential")  diffusion  of  species  or  enthalpy  has  occurred  and  the  elemental  mass 
fractions  and  the  enthalpy  are  linearly  related  to  5  and  to  one  another. 

In  practical  flames  diffusivlties  are  seldom  equal  and  resort  is  made  to  the 
argument  [5]  that  if  the  Reynolds  number  of  the  turbulence  is  high  enough  molecular 
effects  will  be  confined  to  the  high  wavenumber  end  of  the  spectrum.  There  will  be 
no  significant  effect  on  the  main  quantities  of  interest,  such  as  the  swan 
concentrations  and  temperature,  and  the  variances  and  covariances  of  the 
concentration  and  temperature  fluctuations;  these  quantities  are  associated  with  the 
low  wavenumber,  "energy-containing"  range  of  the  turbulence  spectrum.  As  a 
consequence  the  mixtuce  fraction  concept  remains  valid  as  far  as  the  main  features 
of  theoretical  models  are  concerned. 
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In  laboratory  flames  turbulence  Reynold*  numbers  are  seldom  very  high.  In  non¬ 
reacting  Clows  in  the  laboratory,  turbulence  Reynolds  numbers  are  relatively  modest, 
but  they  are  in  general  sufficient  to  ensure  a  separation  of  the  energy-containing 
and  dissipation  ranges  of  the  turbulence  spectrum.  At  flame  temperatures,  however, 
the  kinematic  viscosity  is  increased  tenfold  or  more  and  turbulence  Reynolds  numbers 
are  correspondingly  lowered.  Furthermore,  the  diffusivities  of  some  species  such  as 
molecular  and  atomic  hydrogen  are  so  much  higher  than  the  kinematic  viscosity  that 
it  can  be  expected  that  molecular  transport  will  have  even  greater  effects  on  the 
low  wavenumber  end  of  the  spectrum  for  hydrogen  containing  species.  Thus  one  might 
expect  to  find  significant  affects  of  molecular  transport  on  means  and  variances  of 
quantities  such  as  concentration  and  temperature  in  laboratory  flames. 

Recent  measurements  by  Drake  et  al.  [6,7]  in  a  hydrogen  diffusion  flame  of 
modest  Reynolds  number  show  a  significant  departure  from  what  is  expected  from  equal 
diffusivity  theory.  They  have  used  pulsed  laser  Raman  spectroscopy  to  obtain 
simultaneous  measurements  of  temperature  and  nitrogen  and  hydrogen  concentration 
with  excellent  spatial  and  temporal  resolution.  Correlation  plots  of  concentration 
versus  temperature  for  these  two  species  show  considerable  departure  from  the  curve 
that  is  obtained  assuming  equal  diffusivities,  adiabaticity  and  chemical 
equilibrium.  Predicted  departures  from  adiabaticity  and  chemical  equilibrium  are 
not  qualitatively  or  quantitatively  consistent  with  the  measurements.  The 
measurements  do  in  general  lie  between  the  equal  diffusivity  curve  and  calculations 
made  by  Hiller  and  Ree  [8]  for  a  laminar  flame.  Differential  diffusion  effects  are 
significant  in  the  laminar  flame  and  Drake  et  al.  [6]  consider  that  such  effects  are 
the  explanation  for  their  measurements. 

A  theoretical  medal  for  predicting  differential  diffusion  effects  in  turbulent 
diffusion  flames  is  presented  here.  It  is  used  to  predict  the  departures  from  the 
equal  diffusivity  theory  for  the  conditions  of  the  Drake  et  al.  [6]  flame.  It  could 
also  he  used  to  predict  the  effects  of  differential  diffusion,  or  Reynolds  number, 
on  nitric  oxide  formation  in  such  flames,  a  question  which  has  been  raised  in 
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earlier  work  [9] .  other  uses  of  such  a  model  include  the  prediction  of  the  extent 
of  differential  diffusion  between  the  fuel  and  seed  particles  put  into  the  flow  to 
act  as  a  fuel  tracer  in  the  so-called  marker  nephelometry  technique  [10, 11].  Such 
studies  are  planned  for  the  future. 

In  earlier  papers  [12,13]  the  theory  for  differential  diffusion  in  non-reacting 
flows  has  been  developed.  An  equation  is  derived  for  a  variable  z  which  is  the 
difference  between  mixture  fractions  defined  for  various  species.  Ibis  parameter  is 
then  a  direct  measure  of  differential  diffusion  effects  and  has  zero  values  in  both 
the  unmixed  fluids.  Its  only  source  is  a  source  term  involving  the  difference, 
d,  in  the  diffusion  coefficients  for  the  species  and  the  mean  gradient  of  a  mixture 
fraction.  It  is  shown  that  a  characteristic  scale  for  the  amount  of  differential 
diffusion  Zq  is  related  to  the  characteristic  scale  of  the  mixture  fraction 
and  the  turbulent  diffusivity  dt  by 

*0  !Dt  +  D)  *0  ( 1  * 

where  D  is  the  actual  molecular  diffusivity  for  the  species.  Since  the  turbulent 
diffusivity  increases  with  turbulence  Reynolds  number  Refc  ( =  u'L/v  where  u'  is 
the  rms  velocity  fluctuation,  L  the  integral  length  scale  of  the  turbulence  and 
v  the  kinematic  viscosity)  the  asymptotic  behaviour  at  high  Reynolds  number  is 

zQ  «  1/Ret 

which  is  consistent  with  our  notions  that  the  effects  will  only  be  important  at  low 
Reynolds  number. 

In  chemically  reacting  systems  this  approach  is  complicated  by  the  fact  that 
conserved  scalars,  such  as  elemental  mass  fractions  and  enthalpy,  are  dependent  on 
many  molecular  species  concentrations  and  the  temperature,  each  with  its  own 
effective  diffusivity.  Moreover  the  relative  proportions  of  the  species  change  and 
it  is  not  possible  to  directly  derive  a  simple  weighting  for  these  diffusivlties. 

It  was  thought  earlier  [12]  that  this  problem  could  be  overcome  by  an  appropriate 
scaling  in  the  reaction  zone,  but  this  approach  has  not  proved  fruitful.  In  this 
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paper  a  different  approach  Is  adopted.  Molecular  species  concentrations  and 
teaperature  are  related  directly  to  the  mixture  fraction  (the  equal  diffualvity, 
adiabatic,  chemical  equilibrium  relations)  with  linear  perturbations  for  the  effects 
of  differential  diffusion.  It  is  then  possible  to  define  effective  diffuaivities 
and  differential  diffusivities  for  aleaental  species  and  enthalpy  and  to  formulate 
equations  for  their  differences  or  perturbations  in  their  mixture  fractions. 

In  the  next  section  the  theory  is  dsveloped.  This  is  followed  by  presentation 
of  the  effective  diffusivities  and  other  properties  computed  for  the  hydrogen-air 
system.  The  modeled  turbulence  equations  are  then  solved  numerically  for  the  Drake 
et  al.  (16]  flame  conditions  and  the  results  compared  with  experiment. 

Theory 

In  the  absence  of  significant  body  fores  and  viscous  contributions  to  the 
energy  balance,  the  species  balance  and  energy  equations  may  be  written  [14.15] 

DYi 

p5r“"i-7-id  121 

p  Dt  *  "  V^+  Dt  ’  p*r  (3) 

where  and  w^  are  the  mass  fraction  and  volumetric  reaction  rate  of  species 

1  and  ^  its  mass  flux  relative  to  the  overall  mass  flux  of  the  fluid,  p  is  the 

density,  P  the  fluid  pressure  and  h  the  enthalpy  (including  chemical 

contribution)  given  by 


N  T 

h  =  l  r  (h°  *  J  c  «r)  s  l  rLh 

•  -  rjl  i--«  *  * 
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(4) 
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where  and  h°  are  the  specific  heat  at  constant  pressure  and  the  enthalpy  of 

formation  of  species  i,  and  the  rate  of  radiative  loss  per  unit  mass, 

respectively.  Neglecting  thermal  diffusion  contributions  to  the  heat  flux  vector, 
this  may  be  written 
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*1  “  “  X7T  +  J  ^h 
i-1 


(5) 


with  X  the  thermal  conductivity  and  T  the  temperature.  The  dilute  species 
approximation  is  adopted  for  the  diffuaional  fluxes  ^  so  that  they  have  the 
Pick's  Law  form 


i-1  to  N-1 


N-1 


*  I 

k-1 


‘"k'V 


i  -  N 


(6) 


where  is  the  diffusivity  of  species  i  into  the  diluent  species  N,  here 

assumed  to  be  nitrogen.  Note  that  equations  (6)  preserve  continuity*  this  is  not 
achieved  if  a  Pick's  Law  approximation  is  used  for  the  diluent  species  N.  On 
substituting  equation  (6)  Into  equation  (5)  we  obtain 

N-1 

<1  -  -xvr  +  l  pD  (h  -  hiw.  . 

i»1  1  x  1 

Using  this  and  equations  (6)  in  equations  (2)  and  (3)  and  neglecting  the  substantial 
derivative  of  the  pressure  we  obtain  self  consistent  equations  for  species  and 
enthalpy. 


DYi 

p  —  -  V-lcD^) 
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i  *  N 


(7) 


p  {£  -  7.{X7T  +  l  PDi(h1  -  h^VfJ  -  -p4r 


N-1 


(8) 


The  assumptions  lying  behind  equations  (7)  and  (8)  are  many,  but  they  are 
usually  necessary  to  make  combustion  problems  tractable  [14].  Possibly  the  most 
serious  approximations  in  the  present  case  are  the  use  of  Pick's  Law,  equations  (6), 
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in  a  multicomponent  mixture  and  the  neglect  of  thermal  diffueion.  As  ie  discussed 
later  it  may  be  possible  to  rework  the  problem  with  these  restrictions  removed.  The 
form  of  equation  (8)  is  unusual  but  it  is  the  correct  form  consistent  with  the 
Pick's  taw  assumption  for  dilute  species. 

Element  mass  fractions,  Za,  for  element  m  may  be  defined 


5  l  Ha  iYi 

i-1  a,i  1 


(9) 


where  p  is  the  mass  fraction  of  element  m  in  species  i.  Equations  (7)  may 
m,  l 

be  weighted  by  p  ,  and  summed  to  yield 

St;  i 
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since 
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as  elements  are  conserved  under  chemical  reaction. 

He  consider  a  two  source  (or  stream)  mixing  problem  in  which  the  unmixed 

sources  (or  inlet  streams)  are  denoted  by  the  subscripts  1  and  2.  We  have  then 

Zm(x,t)  -  Zm  2  +  (t  +  V<ZB,1  "  Zm,2)  (11) 

h  -  h2  +  C(h1  -  h2}  +  *h®  l12) 

where  £(x,t)  is  a  mixture  fraction  based  on  one  of  the  elements  (so  that  z_  -  0 

for  that  element)  and  z  (x,t)  and  z.  (x,t)  are  perturbations  for  element  m  and 

m  ~  n  ~ 

the  enthalpy.  In  the  absence  of  differential  diffusion  zm  and  z^  are  zero  and 
E  is  the  universal  mixture  fraction.  In  some  applications  of  interest  h 1  »  h2 
and  so  that  hy  -  h2  is  not  an  appropriate  scale  for  the  enthalpy  perturbation 


z^.  The  scale  Q  may  be  arbitrarily  chosen  and  is  taken  here  as  the  heat  of 
combustion  per  unit  mass  of  stoichiometric  combustion  products. 


The  idea  is  to  substitute  equations  (11)  and  (12)  into  equations  (10)  and  (S) 
respectively  and  derive  equations  for  the  perturbations  ^  and  as  well  as  the 
aquation  for  (.  This  can  be  done  for  a  system  with  any  number  of  elements.  Here 
we  shall  confine  ourselves  to  the  hydrogen/ oxygen/ nitrogen  system.  He  take  hydrogen 
element  as  the  basis  for  £  and  then  only  zQ  for  oxygen  need  be  determined  since 
the  mass  fraction  of  nitrogen  element  is  obtained  by  difference.  The  mixture  is 
then  described  entirely  by  £  and  zQ>  which  determine  its  elemental  composition, 
and  by  z^  which,  together  with  £,  determines  its  enthalpy.  If  the  chemical 
kinetic  rates  are  very  fast  compared  to  the  rate  of  turbulent  mixing  then  the 
molecular  species  and  temperature  will  be  everywhere  the  chemical  equilibrium  values 
for  that  elemental  composition  and  enthalpy.  This  will  be  closely  approximated  for 
the  hydrogen/oxygen/nitrogen  system  of  interest  here  [5,16].  We  have  then  that 

-  Y°(5.*0.«h)  (13) 

T  -  T°(£,z0,*h)  (14) 

where  the  superscript  o  denotes  chemical  equilibrium.  Also 
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These  may  be  substituted  into  equations  (10)  and  (8)  along  with  equations  (11)  and 
(12)  to  yield  after  some  manipulation 

p  H  -  V.tpD^V^)  -  V^pD^Vz^  +  V^pD^Vz^  <17) 

Dzo 

P  — -  V.lpD^)  -  V.(PD0575)  +  ?.<pDohV  (18) 
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°5T-  v'<pDhh7V  “  V’“*V5>  +  7-,pDh0VV  •  p*  (,9> 
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In  the  above 


Vi  “  UH,i/lZH,1  ”  ZH,2’ 

Vi  ■  Vi/(Zo,i  -  zo,2> 

and  the  partial  derivatives  are  evaluated  at  the  local  values  of  C,  *Q  and  zh. 

To  reduce  the  size  of  the  computational  task  we  here  introduce  some 
linearizations.  He  assume  that  zQ  and  are  small  so  that  extraneous  terms  in 

equations  < 17 ) — ( 19 )  may  be  dropped  yielding 


err 


-  V.<oDrrV5)  -  0 


p  5t  “  7*(pDooV  -  7*“*ocvo 


p5T-  7‘<pDhhV  '  7*(pDh?ve)  -  p* 


and  tha  partial  darivativea  appearing  in  equations  (20)  are  then  evaluated  at 
5  «  E,  *«  *  0(  *jj  *  0.  Furthermore  equations  (13)  and  (14)  can  be  linearized 


3Y°  W° 

Hi  -  *°( 5.0,0)  +  aT  *o  +  IT  *h 

0  h 


T  -  T°(  5,0,0)  +  *0  ♦  !r  *h 

0  h 


with  the  derivatives  here  also  evaluated  at  5  ■  6»  “  °»  *j,  “  °* 

Equations  <21)-(23)  are  of  tha  same  fora  as  equations  (6)  and  (7)  of  the 
earlier  paper  (13]  and  may  be  averaged  and  modeled  in  a  similar  way  to  that  paper. 
We  use  density  weighted,  or  Favre,  averaging,  the  axisymaetric  boundary  layer  form 
of  the  equations  and  the  k  -  e  (17]  model  for  the  turbulence.  Equations  for 
5*  s"2,  *Q,  *q2»  *h<  *£2«  5"*q,  5"*J|,  z“z“  result: 


p"  *  +  ^  i  '  r  i?  rpr(-^  ♦  ®K)  ■  0 


*  +  ‘pv~  ^  -  7  I?  +  ”55 }  ^  "  _pCWf )  -  Cg2p  k  5*2 
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with  y  being  the  ratio  of  the  equilibriua  specific  heats*  It  is  found  that  bv 
is  largely  independent  of  £.  The  Pavre  averaging  of  vs  and  ay  is  done  in  the 
normal  manner  using  the  Pavre  probability  density  function  for  £. 

Boundary  conditions  for  z^,  z^2,  £"z“ ,  and  zgz£  are  that  they  have  zero 
value  in  the  two  inlet  streams  and  a  zero  normal  derivative  at  the  axis  of  symmetry. 
Non-zero  values  for  zfc  arise  from  the  source  term  on  the  right  hand  side  of 
equation  (27)  and  gradients  in  z^  give  rise  to  the  variances  and  covariances. 
Property  Data  for  the  Hydrogen/Alr  System 
The  diffusion  coefficients  appearing  in  equations  (21)  to  (23)  and  defined  in 
equations  (20)f  were  determined  in  the  following  manner.  The  species  diffusivities 
were  assumed  to  be  those  into  nitrogen  and  these  were  obtained  from  power  law  fits 
to  the  Chapman-Knskog  theory  derived  by  Mitchell  [19] 

Di  *  Di,N2  ”  fT1,67/P  cm2/s  . 

Similar  fits  to  the  thermal  conductivity  of  component  gases  as  derived  by  Mitchell 
[19]  were  also  used: 

”  cTa  cal/cm-s-K  . 

These  were  mole  fraction  weighted  to  obtain  the  overall  thermal  conductivity. 
Molecular  viscosities  were  also  determined  for  each  species  from  Mitchell' s  fits 

Ul  -  oT0,6756  g/cm-s 

and  the  overall  viscosity  (used  in  the  momentum  equation  model)  determined  by  mole 
fraction  weighting.  In  all  cases  T  Is  in  Kelvin  and  the  pressure  P  is  in 
atmospheres.  The  constants  for  the  fits  are  given  in  Table  1.  The  mole  weighting 
of  viscosity  and  thermal  conductivity  are  gross  approximations  but  are  consistent 
with  the  treatment  of  the  mass  diffusivities. 

The  density,  temperature  and  composition  were  obtained  as  a  function  of  the 
mixture  fraction  (mass  fraction  of  hydrogen  fuel.  *0  and  r..n  zero)  by  calculation 
for  adiabatic,  equilibrium  reaction  with  dry  air  using  the  program  of  Gordon  and 
McBride  [20]  and  the  derivatives  with  respect  to  £  being  determined  by  finite 
difference  techniques.  Enthalpies  of  individual  species  were  obtained  from  the 


Table  It  Conatanta  used  In  fita  for  diffusivlty,  thermal 
conductivity  and  vlacoalty  [19] 


fT1,67/P, 

cm2 /si  lt  - 

cTa,  cal/cm-a-Ki  ^  ■ 

er°-«756,  ,/, 

Species  i 

f  x  105 

c  x  107 

d 

e  x  10* 

H 

8.171 

137.2 

0.6484 

1.5426 

«2 

5.525 

51.84 

0.7681 

1.8405 

OH 

2.030 

14.804 

0.7601 

4.115 

h2o 

1.905 

0.8304 

1.1748 

0.5083 

**2 

1.488 

7.6893 

0.7722 

3.6974 

O 

2.076 

22.219 

0.6547 

4.0387 

°2 

1. 523 

7. 1352 

0.7968 

4.4203 

JANA?  tables  [21].  Partial  derivatives  with  respect  to  zQ  and  a,,  required  for 
the  derivatives  in  equations  (20),  (24),  (25)  and  (34)  and  (35)  were  obtained  by 
finite  difference  techniques  after  running  the  Gordon  and  McBride  [20]  adiabatic 
equilibrium  calculation  for  several  values  of  the  oxygen  content  and  the  temperature 
of  the  air.  In  one  case  this  was  not  necessary  as 


(36) 


where  c*  la  the  equilibrium  specific  heat  of  the  mixture.  The  enthalpy  scale 
factor  Q  is  arbitrarily  chosen  as  the  heat  of  reaction  per  unit  mass  of 
stoichiometric  mixture  and  is  equal  to  811.4  cal/g. 

The  diffuaivities  D„„,  D„,  and  are  shown  in  Figure  1.  Algebraic  fits 

to  these  as  used  in  the  numerical  integration  are  shown  in  Table  2.  It  can  be  seen 
that  Df^,  which  can  be  Interpreted  as  the  effective  diffusivlty  for  hydrogen 
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Figure  1.  Effective  diffusivities  for  hydrogen  and  oxygen  elements. 


Table  2i  Algebraic  fits  to  property  data 


Property 

Units 

Range  of  £ 

Formula 

D« 

o»2/a 

0  <  £  <  0.0325 

0.0325  <  £  «  1 

-  0.26  +  257 £  +  3.5x108£5 

-  3.8  -  3.03£  ♦  27.5(1.0001  -  £)10 

°oe 

CM2/* 

0  <  £  <  0.025 

.025  <  £  <  1 

DQ£  -  -1.83  -  2000 £ 

D0?  -  -0.51  -  0.385(£-1'3  -  1) 

°oo 

c*2/e 

0  <  £  <  0.03 

Dqo  -  0.21  +  206. 7  £  1778£2 

0.03  <  £  <  1 

Oqq  -  1.261  -  £  +  9.215(1.0001  -  £)9*5 

V 

c»2/a 

0  <  £  <  0.020 

11^-  -2.5  -  575£  -  52,800 E2  ♦  1.0495x1010 £5 

0.02  <  £  <  0.035 

0.035  <  £  <  1 

Dh?  -13  +  145  exp {-14860 ( £  -  0.033)2} 

Dh?-  -3.93  -  58(1.0001  -  £) 3  +  1.25£-1*475 

Dhh 

cm2/ a 

0  <  £  <  0.048 

Dj^  -  0.214  ♦  150£  ♦  8182 £2  -  5.096x10?£5 

0.048  <  £  <  1 

0^  -  7.39  -  6£  ♦  9.88(1.0001  -  £>5'34 

V 

c«2/e 

0  <  £  <  .0265 

V  -  0.154  ♦  150. 9£ 

.0265  <  £  <  1 

V-  1.96  -  0.9£  ♦  2.68(1.0001  -  £)6*8 

v. 

*3/kg 

0  <  £  <  0.0275 

-  0.845  +  335. 3£  -  2830£2 

0.0275  <  £  <  0.4 

va  -  6  +  6.36( £  -  0.0245)0'208 

0.4  <  £  <  1 

vm  -  10.4  +  1.76£ 

*v 

«3/kg 

0  <  £  <  0.03 

Sy  -  0.0254  +  4£  1.925x107£5  -  6.173x1021 

0.03  <  £  <  1 

Sy  -  -9.8  +  0. 1314( £  +  0.025)"1*39 

♦ 

a’1 

0  <  £  <  0.033 

♦  -  2.97  exp{-6172.8(0.031  -  £)2) 

0.033  <  £  <  1 

♦  -  12.3(1  -  £2)6,7 
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element,  reaches  quits  high  vsluss  just  on  the  rich  side  of  stoichiostetric  which 

occurs  at  E  «  E  »  0.0283.  It  is  seen  that  on  the  lean  side  of  stoichiometric  this 
a 

diffusivity  is  2  to  5  times  the  kinematic  viscosity  (shown  in  Figure  3  and  Table 
2).  This  low  effective  molecular  Schmidt  number  may  have  implications  for  the 
turbulence  structure  in  the  outer  part  of  the  flame  at  low  Reynolds  number.  The 
differential  diffusivity  can  be  interpreted  as  the  difference  between  the 

diffusivity  of  oxygen  and  hydrogen  element.  It  is  negative  over  the  whole  range, 
indicating  that  the  diffusivity  of  oxygen  element  is  lower  than  that  for  hydrogen 
element.  More  than  that,  it  can  be  seen  that  on  the  lean  aide  of  stoichiometric 
(£  <  0.033,  in  fact)  the  effective  diffusivity  for  oxygen  element  (dq^  +  D , is 
negative,  indicating  that  oxygen  element  diffusion  is  counter  to  the  gradient.  This 
arises  from  the  fact  that  the  diffusivity  of  HjO  is  greater  than  that  for  Oj. 

The  radicals  H  and  OH  make  very  significant  but  not  dominant  contributions  near 
stoichiometric.  The  neglect  of  superequilibrium  radical  concentrations  (due  to 
finite  rate  kinetics  [16])  may  be  a  significant  overall  contribution  at  moderate  to 
low  Demkohler  number. 

The  diffusivities  and  D^h  are  shown  in  Figure  2  and  algebraic  fits  to 

them  are  given  in  Table  2.  It  is  seen  that  D.  ,  is  positive  near  stoichiometric 

ht 

but  is  negative  outside  this  range.  This  complex  behaviour  should  give  rise  to  some 
interesting  effects  on  the  flame  temperature.  The  values  of  the  diffusivltiea  shown 
are  very  large  compared  with  the  kinematic  viscosity.  It  should  be  remembered  that 
the  scaling  for  and  ia  arbitrary  since  it  depends  on  the  value  of  Q. 

However,  the  choice  of  Q  as  the  heat  release  per  unit  mass  of  combustion  products 
should  mean  that  these  diffusivities  are  the  effective  thermal  diffusivity  operating 
on  the  flame  to  air  temperature  difference. 

xgure  3  shows  the  kinematic  viscosity,  v,  and  the  specific  volume 
functions  vg  and  ay  defined  in  equations  (32)  and  (34).  Algebraic  fits  are  once 
again  shown  in  Table  2.  It  is  seen  that  ay  undergoes  a  very  strong  variation,  its 
large  negative  value  on  the  rich  side  of  stoichiometric  being  due  to  the  very  strong 
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Kinematic  viscosity,  specific  volume  and  influence  coefficient 
of  oxygen  differential  diffusion  on  specific  volume. 


•f fact  of  a  decrease  of  elemental  oxygen  on  tha  flams  temperature.  On  tha  othar 
hand  by  was  found  to  vary  froa  9.57  m3/kg  at  5*0  and  1,  dipping  to  a 
■tin in ua  of  5.1  mVkg  at  5  “  0.03.  Such  variation  la  of  second  order  and  an 
average  value  of  by  ■  7.3  a3/ kg  has  been  used  in  the  calculation. 

The  radiative  loss  rate  term  6  was  calculated  froa  the  absorptivity  of  water 
vapour  [22]  using  X°  and  T°  assuming  the  flaae  to  be  optically  thin,  h 

HjO 

functional  fit  with  5  has  been  determined  as  shown  in  Table  2. 

In  converting  the  computational  results  in  teras  of  cQ  and  to  the 

effects  on  species  concentration  and  temperature ,  it  is  necessary  to  use  the 
functions  and  derivatives  that  appear  in  equations  (24)  and  (25).  Figure  4  gives 
these  for  temperature  and  Figures  5  and  6  for  nitrogen  and  hydrogen  respectively. 
Results  for  other  species  can  also  be  obtained  as  well  as  for  combined  variables 
such  as  nitric  oxide  formation  rate. 

Results 

numerical  solution  of  equations  (26)  to  (28)  have  been  carried  out  using  the 
finite  difference  computer  program  and  methodology  reported  earlier  [18]. 

Turbulence  modal  constants  have  been  kept  the  same  and  the  turbulent  Schmidt  number 
is  taken  as  0.7.  Tha  special  model  for  low  Reynolds  number  turbulence  of 
Jones  and  Launder  [17]  has  been  included  in  the  modeled  equations  for  velocity, 
turbulence  kinetic  energy  and  dissipation  which  are  solved  simultaneously.  The 
solution  procedure  uses  a  simple  model  [18]  for  the  probability  density  function  of 
(  involving  a  clipped  gaussian  distribution  plus  an  intermittent  spike.  The 
intermittency  is  found  from  an  empirical  correlation  dependent  on  the  mean  and 
variance  of  £.  This  pdf  is  used  to  weight  the  property  data  shown  in  Table  2  to 
obtain  the  necessary  averages  D^,  DQ^,  v,  p  etc. 

Calculations  have  been  made  for  the  flame  of  Drake  et  al.  [6]  in  which  a  jet  of 
hydrogen  3.2  mi  diameter  exits  with  a  mean  velocity  of  50  a/ s  into  a  co- flowing 
stream  of  air  with  a  velocity  of  10  m/».  No  information  is  available  on  the  initial 
velocity  or  turbulence  profiles  or  the  axial  pressure  gradient  in  the  tunnel.  Since 
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the  Reynolds  number  is  low  at  1500  it  has  baen  assumed  that  tha  hydrogen  flow  at  jat 
ax it  has  a  parabolic  valocity  profila.  Tha  axtarnal  velocity  and  turbulence 
intanaity  profila  has  baan  arbitrarily  ohosan  as  a  scalad  varsion  of  that  measured 

by  Rambach  at  al.  [23]  in  a  similar  apparatus*  Tha  dissipation  langth  seals 
3/2 

l*e  s  k  /c  is  arbitrarily  chosan  at  0*2  tiaas  tha  jat  diaaatar.  Turbulanca 
intansitias  inside  tha  jat  are  also  sealed  frost  Rambach  at  al*  [ 23]  with  a 
dissipation  langth  seals  of  0.1  tiaas  tha  nossla  disaster.  Tha  axial  pressure 
gradient  ia  chosan  at  -8  Fa/a  which  causes  a  2*  acceleration  of  the  axtarnal  flow  at 
100  nossla  disasters  downs traaa.  Computation*  have  baan  aada  with  and  without 
radiative  loss* 

Figure  7  shows  results  for  %,  s^,  and  sQ  along  the  centreline  of  this 

flaas*  It  is  seen  that  cQ  and  s^  are  quite  large  and  that  tha  linearisations  of 

aquations  (21)-(26)  and  (33)  are  not  justified  in  this  flaas.  (Such  linearisations 

will  of  course  be  justified  in  flaaas  of  higher  Reynolds  nuaber  where  sQ  and  s^ 

will  be  auch  saallsr.)  It  is  seen  that  s^  is  positive  near  tha  jet  exit  and  is 

negative  further  down* tr earn.  This  can  be  ascribed  to  the  fact  that  ,  near  the 

hi 

centreline  is  negative  near  the  jet  and  positive  downstream.  Figure  8  shows  the 

radial  variation  at  x/b  of  55.  (o  is  here  jot  diameter  and  a  -  D/2.)  Table  3 

lists  the  various  diffusivitias  across  the  flow  at  x/b  •  55. 

The  change  of  sign  of  zQ  in  the  outer  part  of  the  profile  is  related  to  the 
2~  2 

change  in  sign  of  )  (/9r  and  is  necessary  to  conserve  a  zero  total  flux  of  zQ 
when  integrated  across  the  flow.  This  necessary  condition  forma  a  good  check  on  the 
numerical  accuracy  of  the  computation  and  in  our  case  the  flux  of  zQ  was  less  than 
10"3  of  the  flux  of  5  up  to  x/D  of  100.  A  similar  result  for  this  check  is 
obtained  for  z^  in  the  no  radiation  case.  The  estimation  formula,  equation  (1), 
with  d  here  as  D^,  0  ”  and  -  \|j/0.7  gives  order  of  magnitude  estimates 

for  zQ  and  z^  on  the  centreline  with  a  constant  of  proportionality  of  0.5  and 
0.2  when  the  diffusivitias  are  calculated  in  the  middle  of  the  shear  layer 
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Figure  7.  Results  for  mixture  fraction  and  its  perturbations  on 
centerline  of  hydrogen/air  diffusion  flame  of 
Drake  et  al .  [6] . 
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(€/{q  “  0.5)  and  ^  is  tha  valua  of  £  at  tha  cantralina.  Thaaa  can  ba 
compared  with  tha  valua  of  0.5  found  in  tha  much  alaplar  nonreacting  flow  casa  [13]. 


Tabla  3i  Naan  diffuaivitiaa  and  othar  var labia#  aoroaa  tha 
flame  shaar  layar  at  x/D  -  55 


r/a 

0 

2.95 

4.44 

6.27 

8.13 

1.0 

0.788 

0.581 

0.279 

0.028 

(u  -  U#)/(UC  -  u#) 

1.0 

0.584 

0.118 

-0.052 

-0.044 

«2/. 

15.4 

17.8 

20.2 

9.3 

0.86 

Doe  cm2/. 

-10.5 

-14.5 

-21.5 

-41.8 

-6.5 

°oo  c-2/* 

5.3 

6.1 

7.0 

5.9 

0.71 

Dhe  cm1/. 

4.1 

24.8 

63.5 

32.5 

-4.2 

^hh 

13.2 

13.9 

14.0 

7.6 

0.63 

v  cm2/. 

3.4 

3.6 

3.9 

3.5 

0.51 

vT  cm2/. 

8.4 

3.9 

0.18 

0.02 

0.03 

Rat 

67.8 

49.7 

14.3 

1.89 

7.1 

Comparison  of  tha  computed  flame  flald  with  tha  measurements  of  Oraka  at  al. 

[6]  indicate#  that  the  computed  flame  is  about  ten  percent  longer  than  the  actual 

flame.  This  discrepancy  may  be  due  to  the  mismatch  of  initial  conditions,  assumed 

versus  actual.  The  agreement  obtained  is  really  quite  good  considering  that  tha  low 

Reynolds  number  affects  on  the  turbulence  are  quite  large.  Table  3  shows  that  in 

2 

the  middle  of  the  shear  layer  the  turbulence  Reynolds  number  (Ret  =  k  /ve)  is 
about  40.  This  reduces  C^,  the  coefficient  in  the  turbulent  viscosity  model  by  a 


-25- 


factor  of  about  3  with  consequent  reduction  of  turbulence  levels  and  the  rate  of 
f lane  spread. 

The  effects  of  differential  diffusion  on  the  molecular  species  concentrations 
and  temperature  could  be  obtained  by  weighting  the  linearized  equations  (23)  and 
(24)  by  the  joint  pdf  for  £,  zQ  and  s^.  This  would  give  average  temperatures  and 
concentrations  and  average  departures  from  the  equal  diffusivity  relations 
represented  by  the  lead  terms  in  those  equations.  The  computations  carried  out  to 
date  do  not  give  us  enough  information  to  compute  this  joint  pdf.  For  our  present 
purposes  it  will  be  sufficient  to  calculate  characteristic  temperature  and  mole 
fraction  departures  defined  by 


(37) 


(38) 


the  influence  coefficients  3T°/  JZj^  and  aK°/3zfc  being  computed  at  £•  Figures  9 
and  10  show  these  characteristic  temperature  and  mole  fraction  departures  on  the 
centreline  and  across  the  flame  at  xA>  >40.  It  is  seen  that  the  temperature 
departures  are  very  large  indeed.  They  are  an  order  of  magnitude  larger  than  the 
average  temperature  departures  found  for  chemical  non-equilibrium  [16].  It  is 
apparent  that  the  departures  will  still  be  substantial  in  the  flame  that  we  have 
been  studying  at  Sydney  Uhiversity  [24,11,9].  There  could  be  significant  effects  on 
nitric  oxide  formation.  Figure  10  shows  the  breakdown  of  the  temperature  departure 
into  its  components  arising  from  the  differential  diffusion  of  oxygen  element 
dCT0,  the  first  term  in  equation  (37)  and  that  due  to  enthalpy  ACT^,  the  second 
term  in  equation  (37).  It  is  seen  that  oxygen  element  diffusion  is  the  dominant 
feature  on  the  rich  side  of  the  flame.  For  the  species  mole  fractions  the  effect  is 
almost  all  due  to  the  oxygen  element  differential  diffusion. 
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Figure  9.  Characteristic  departures  of  temperature  and 
molecular  nitrogen  and  hydrogen  mole  fraction 
due  to  differential  diffusion  on  flame 
centerline. 


the  flame  at  x/D 


Figure*  11  *nd  12  show  th***  result*  plotted  as  correlation  diagrams  between 
temperature  and  concentration.  Th*  differential  diffusion  theory  curves  shown  are 
for  temperature  computed  as 

T(«>  i  T°(5)  ♦  A°*(C>  (39) 
cross  plotted  versus  concentration  computed  as 


ctU) 


n  «  Q  v 

9  *llV  * 

~  *  O'*  o  .* 

R  *  <C)  ♦  i  t(() 


(40) 


where  9  is  the  pressure  and  R  the  universal  gas  constant.  Th*  equal  diffusivity 
theory  curves  are  then  simply  equations  (39)  and  (40)  cross  plotted  with  ACT  and 
ACX^  set  equal  to  sero.  Th*  cross  plotting  is  done  from  th*  data  obtained  across 
the  flasM  at  x/D  -  SS.  This  is  shown  compared  with  th*  deta  of  Drake  at  el.  [0] 
at  x/D  ■  50  which  is  th*  comparable  station  ainoe  the  computations  give  a  longer 
flam*.  The  cross  plot  correlation  is  not  very  sensitive  to  x/D  in  this  part  of 
th*  flams.  It  can  be  seen  that  th*  agreement  i*  in  general  very  good,  th* 
differential  diffusion  theory  lying  much  closer  to  th*  middle  of  th*  swasured  data 
than  for  the  equal  diffusivity  theory.  The  wide  scatter  of  th*  data  is  due  partly 
to  the  Boisson  statistics  of  th*  measurement  method  but  also  due  to  turbulence.  In 
th*  absence  of  differential  diffusion  and  measurement  errors  all  points  should  lie 
on  the  equal  diffusivity  curve.  With  differential  diffusion  such  perfect 
correlation  would  only  occur  if  th*  fluctuations  in  sQ  and  s^  are  perfectly 
correlated  with  %.  It  is  proposed  to  calculate  these  correlation*  and  compare  them 
with  the  more  recent  measurements  of  Drake  et  al.  [7}« 

Although  the  agreement  between  theory  and  experiment  is  most  encouraging  this 
should  be  viewed  against  improvements  needed  in  the  theory  to  obtain  really 
satisfactory  predictions  in  this  case.  The  linearizations  involved  in  equations 
(21)— (23)  could  easily  be  eliminated  but  those  involved  in  equations  (24)  and  (25) 
and  the  effective  dif fusivitles  are  not  so  easily  handled.  Th*  Tick’s  Law 
assumption  could  be  eliminated  by  solving  the  Stefan-Maxwell  equations  by  using 
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Figure  11.  Predicted  effect  of  differential  diffusion  on  correlation  < 
temperature  and  hydrogen  concentration.  Experimental  data 
of  Drake  et  al.  [6] . 
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Figure  12.  Predicted  effect  of  differential  diffusion  on  correlation  of  temperature 
hydrogen  concentration.  Experimental  data  of  Drake  £t  al.  [6) . 


aquations  (IS),  (  IS)  and  (14)  and  (13)  ao  that  nw  expressions  for  tha  effective 
diffusivltiaa,  aquations  (20) ,  ara  arrivad  at. 

Conclusions 

A  thaory  for  dlffarantlal  molecular  transport  In  turbulent  diffusion  flames  has 
baan  developed  from  fundsmental  principles.  Effective  dlffarantlal  diffusivities 
for  elemental  apacias  and  anthalpias  hava  baan  derived  and  valuas  oomputed  for  tha 
hydropan/air  system.  Thasa  valuas  ara  found  to  ba  auch  la r par  than  tha  klnaaatio 
vlaoosity  which  la  already  large  In  auch  flames.  Solutions  of  tha  turbulant  fora  of 
tha  aquations  for  tha  perturbations  in  alaaant  and  enthalpy  mixture  fractions  yield 
results  which  ara  In  pood  agreement  with  laser  Raman  measurements  In  a  low  Reynolds 
number  flame.  Tha  thaory  could  ba  improved  by  incorporating  low  Reynolds  number 
affects  In  th  turbulence  modal,  eliminating  soaw  of  tha  linearisations  and  avoiding 
tha  Pick's  Law  assumption. 
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